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Abstract 



^ I We present a method for constructing the complete set of meson-exchange 

corrections to baryon observables described by covariant chiral quark models. 
■ The meson corrections are expressed in terms of an unperturbed valence- 

. quark Green's function. This method is illustrated for the case of an NJL 

model. It also enables us to discuss, in terms of approximations used for this 
unperturbed Green's function, the treatments of meson corrections in two 
' other chiral quark models. 

iu: 

^ ■ I. INTRODUCTION 

> : 

^ . Chiral symmetry is one of the most important features of QCD at low energies. This 

^ I approximate symmetry is spontaneously broken by the quark condensate, leading to con- 
stituent masses for the quarks and to pions which are, to a good approximation, the corre- 
sponding Goldstone bosons. As a result pions are much lighter than all other mesons and 
so the cloud of virtual pions surrounding a baryon can make significant contributions to the 
properties of that baryon. Any realistic model for the structure of baryons should therefore 
take account of their pion clouds. 

As approximate Goldstone bosons, pions with low momenta interact relatively weakly 
with other hadrons. This forms the basis for chiral perturbation theory (ChPT), in which 
hadron properties and scattering amplitudes are expanded in powers of momenta and the 
pion mass (for a review, see fl]]). It also provides a motivation for treating the pion cloud 
perturbatively in chiral quark models. Such a treatment ensures that at least the longest- 
range effects of the pion cloud are described correctly, since these involve virtual pions with 
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low momenta. These effects can be checked by comparing contributions to hadron properties 
that depend nonanalytically on the square of the pion mass with the model-independent 
predictions of ChPT. 

A range of chiral quark models for baryon structure can be found in the literature. These 
include chiral or "cloudy" bag and soliton models 0,^] and chiral versions of constituent 
quark models In all of these, pions are introduced as elementary degrees of freedom 

and in many cases their contributions are calculated perturbatively. 

In contrast, models of the Nambu-Jona-Lasinio (NJL) type are based on relativistic 
quarks only, and the quark condensate and pions are generated by the interactions between 
these quarks (for a review, see: As such they represent a step closer to QCD, although 
problems with implementing confinement remain to be solved satisfactorily. In these models 
baryon states have been constructed by solving the relativistic Faddeev equation p|-p!T|, but 
this approach does not include the pion cloud.Q The part of the cloud corresponding to the 
pairwise meson exchange between quarks has been studied in Ref. Recently a covariant 
version of the treatment used in the cloudy-bag model has been applied to calculate pion 
cloud effects based on a Faddeev approach to the nucleon |T6| , p^ . Also, a quite different. 



noncovariant approach to the calculation of pion cloud effects in models with composite 
pions has been proposed by Bicudo et al. JTB| . 

In the present work we address the problem of mesonic corrections within an NJL model, 
deriving covariant expressions for the complete set of one-loop corrections to baryon masses 
and couplings to currents. The approach is based on the "gauging of equations" method 
||T9| for constructing photon couplings to strongly interacting systems of particles. Our 
approach is applicable to any set of one-boson contributions to such systems and was recently 
used to construct a complete set of electromagnetic corrections [^|. As has been noted 



previously in the context of NJL models for mesons pT|-p^ , such completeness is important 
since it ensures that physical quantities calculated in this approach satisfy the relevant 
symmetry constraints. The general features of the approach also allow us to comment on 
the approximations used in cloudy-bag and chiral quark models, which have recently been 
the subject of some debate 



II. COVARIANT MESONIC CORRECTIONS 

A method has recently been presented for calculating all possible electromagnetic cor- 
rections, at order e^, to any quark or hadronic model whose strong interactions can be 



represented nonperturbatively by integral equations [20|. We will use this method to cal 



culate in a covariant manner all mesonic corrections to a relativistic quark model. For this 



^We should mention that a completely different approach to baryons in these models uses a mean- 
field treatment where baryons are described as Skyrmion-like solitons |l2| , |l^ . This is expected to 
be valid in the limit of a large number of colours. However in the realistic case of three colours the 
Faddeev approach is likely to be more appropriate. An attempt to hybridise the two approaches 



can be found in Ref. |14| 
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purpose all we need to do is to summarise the main model-independent results of Ref. pO| , 
applied to exchange of mesons rather than photons. 

In such a model, the Green's function G describing a bare system of quarks satisfies an 
integral equation whose symbolic form is 

G = Go + GqKG, (1) 

where Gq is the "free" Green's function. Note that this need not describe the propagation 
of three free quarks. In a Faddeev approach, for instance, it is the Green's function for one 
free quark and an interacting pair of quarks (often referred to as a diquark, although the 
pair need not be bound). 

The lowest-order mesonic correction to G is denoted by 5G. This should contain all 
possible insertions of the meson propagator, D^^. In chiral quark models with elementary 
mesons this is just a standard boson propagator, but in NJL-type models it is essentially 
the quark-antiquark T-matrix, as discussed below. A four-vector-like notation is used for 
the indices and v which label the meson channels. This will allow for applications to 
models with vector or axial mesons. The propagator is diagonal in its indices, except for the 
mixing between pseudoscalar and axial channels which occurs in models with spin-1 mesons. 
The set of corrections includes mesonic insertions within the kernel K which do not start or 
finish on an external quark line. These provide a contribution 5K to the kernel. Similarly 
insertions within the free Green's function Gq provide a contribution 5Gq. 

All contributions with a meson which connects two elements of the integral equation 
Eq. (|l]) are constructed by the "gauging of equations" method which in this case attaches 
an external meson in all possible ways to each of the unperturbed quantities Gq, K and G. 
The corresponding vertices are given by 

r^ = G^^G^Go\ r'^ = r(; + ir'^, G'^ = Gr'^G. (2) 

An equation for the complete set of lowest-order mesonic corrections to G, denoted by 5G, 
then follows from Eq. (0): 

5G = 6Go + 6G0KG + G06KG + G0K6G + [G^K^G + Gl^KG'' + GqK^'G'') D^,^. (3) 

Writing 

6G = GAG, (4) 

this can be formally solved to give 

A = 6K + G^'SGoGo' + {T^GT" - T'o'CoE^) D^,. (5) 

The quantity A consists of the complete set of one-loop mesonic corrections to the un- 
perturbed kernel K. The third term on the right hand side (RHS) contains the physical 
intermediate states of three quarks plus a meson. There are significant cancellations between 
the second and fourth terms, leaving contributions which correct the Go piece of the third 
term for omissions and overcounting. Note that the contributions of these terms include dis- 
connected pieces where the particles in Go are individually dressed by meson loops. These 
can be easily recast as corrections to the "free" propagator Go. 
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Although our main focus here is on covariant approaches, it is worth noting that the same 
formahsm also applies to time-ordered perturbation theory. The only difference is that the 
Green's functions involved are functions of three-momenta rather than four-momenta. A 
method for applying time-ordered perturbation theory to mesonic corrections in models of 
this type has recently been outlined in Ref. |T^, although that work proposed truncating 
the sum over intermediate baryon states to just the nucleon and A. However, within a 
time-ordered approach, it is difficult to ensure that vertex functions satisfy the constraints 
of Lorentz invariance. This is necessary if one wants to calculate processes involving finite 
momentum transfer to a baryon, or to treat correctly recoil of the intermediate baryons 
in meson-loop diagrams. In addition, a time-ordered approach would make it difficult to 
include quark "z-diagrams" , since these would go beyond a basis of states with three valence 
quarks (plus mesons). Such effects are expected to be important in models where the quarks 
have strong couplings to scalar mesons and couplings to pions with a pseudoscalar form. It 
is also difficult to regularise time-ordered approaches to NJL models in a way that respects 
Lorentz invariance. For these reasons we have concentrated here on developing a covariant 
treatment. 



III. NJL MODEL FOR THE NUCLEON 

The NJL model P] is a relativistic quark model which reflects many aspects of the 
spontaneously broken chiral symmetry of QCD. It therefore provides a useful framework 
within which to investigate the role of mesonic corrections. 

The simplest version of the model is defined in terms of an isospin-doublet colour-triplet 
quark field ip and has the Lagrangian density 

C = ip {i^ - mo) 'ip + [{tptpY - {4^15^^^] (6) 

where r is the vector of isospin Pauli matrices, and mo is the current quark mass (which 
explicitly breaks chiral symmetry). 

Chiral symmetry is spontaneously broken in this model by condensation of quark- 
antiquark pairs. These generate a constituent mass for the quarks, which can be found 
by solving the Schwinger-Dyson (SD) equation for the quark propagator 

S{p) = So{p) + So{pMp)S{p) (7) 

where Sq{p) is the bare quark propagator and the self-energy E is, in the Hartree approxi- 
mation, 

m = -W.J (04tr[5(fc)]. (8) 

Here "tr" denotes a trace over Dirac, colour and flavour indices. The integral over four- 
momenta is, of course, divergent and so needs to be regularised. This may be done either 
by using one of the various cut-off schemes which can be found in the literature, or by 
introducing form-factors to smear out the local interaction. Provided that the coupling g^^ 
is large enough, this SDE has a nonzero solution for S. 
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The quark-antiquark T-matrix can be found by solving a Bethe-Salpeter (BS) equation, 
which in the random-phase approximation is a simple algebraic equation. For example, in 
the pion channel this yields 

where the basic quark-antiquark loop integral is 

TlAp')S^J = J ^tr [z^,nS{p + k)t^,TjS{k)] . (10) 

In the simplest NJL model there is a similar T-matrix in the scalar isoscalar (a-meson) 
channel. In more general versions of the model, interactions in the vector and axial-vector 
channels are treated in the same way. 

In the treatment of the model considered here, the nucleon is described by a three-quark 
vertex function which satisfies a four-dimensional Faddeev equation. The local nature of 
the quark-quark interaction means that the structure of this integral equation is similar to 
that of a BS equation. At lowest order the kernel K can be pictured as exchange of a quark 
between a quark and two interacting quarks (a "diquark"). 

The coupling strengths of the quark-quark interaction can be found by Fierz-transforming 



it into these channels |]I0| . The colour-3 diquark channels are the relevant ones for formation 
of a nucleon. For a local interaction, there are two such channels: scalar isoscalar and axial- 
vector isovector. The transformed interaction Lagrangian has the form 



a. i,a 

(11) 



where 




{n^k = ^h;^a^k, (12) 



and Qs and Qa are the couplings in these channels.^ For the original NJL model, these 
couplings can be related to as described in Ref. In the scalar-diquark channel, the 
propagation of two interacting quarks can be described by the T-matrix, 

where Ils{p) has the same form as in the pion channel, Eq. (PH]). The interaction between 
two quarks in the axial diquark channel leads to a T-matrix whose form is similar to that 
in a vector-meson channel (see Ref. ll^ for more details). 



^ Other contributions including colour-6 diquark channels and the chiral partner of the axial di- 
quark, V'7/i75C'T2/?"V'^, are generated by the Fierz-transformation. We have not written these 
down since they are not needed for the present discussion. 
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FIG. 1. The Faddeev equation for the quark-diquark Green function. 



(a) 





FIG. 2. The mesonic corrections, A, to the nucleon in the NJL model. Wavy hues denote 
meson propagators in the RPA and double lines the corresponding diquark propagators. In line 
(a) the first diagram is the term 5K in Eq. (^), and the rest form Gq^5GqGq^ . The diagrams of 
line (b), along with the ones formed by swapping initial and final states in the last three diagrams, 
make up T^GV^ Dfj^^. Those of line (c) form the subtraction term — FQCoFQ^^jy. In addition 
tadpole-type insertions on the quark lines need to be included. For reasons of space, we show only 
one representative diagram of this class as (d) . 

The corresponding Faddeev equation is represented diagrammatically in Fig. It can 
be written in the form of Eq. (0) by defining the "free" quark-diquark propagator 

G,{P) = S{P-p)D,{p). (14) 

The quark-diquark interaction kernel is given by the quark exchange term 

^a'ab',p)=75/5'^^(p'+p)75/5'^'. (15) 

When the axial diquark channel is included, the kernel must be extended to a 2 x 2 matrix 



which couples the two spin-isospin channels, as described in Ref. |10 . 



IV. MESONIC CORRECTIONS TO THE NUCLEON 

The mesonic corrections to this model can be found by applying the general formulation 
of Eqs. (0)-(|]) to the Faddeev equation outlined above. The resulting expression for A is 
represented diagrammatically in Fig. |^. Mesonic corrections to the nucleon in the NJL model 
have previously been studied by Ishii ||T3[ but it should be noted that the set of diagrams in 
that work is incomplete. Only the diagrams of Fig. ^ were included, with subtracted meson 
propagators as discussed below. 

The full set of diagrams includes not only the meson exchanges between the quarks 
treated by Ishii [|T3| but also the meson clouds of the quark and diquark. In calculating 
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FIG. 3. The parts of the mesonic corrections where the subtracted meson propagator of 
Eq. (Ill) should be used. 



the pion-cloud contributions it is essential to attach the pion in all possible ways to the 
dressed quark propagator obtained from the Hartree approximation ||2l]j2^ . As well as direct 
dressings of the quark lines, as in Fig. |^(a), this includes attachments within the Hartree 
self-energy and generates tadpole diagrams like Fig. H(d) where a pion loop is connected to 
a quark line via an interacting quark- ant iquark pair in the scalar, isoscalar channel. The 
exchanged quark-antiquark pair is represented in Fig. 0(d) by a sigma-meson propagator. 
There is a strong cancellation between these diagrams and the basic pion loop attached to a 
single quark. This cancellation is needed to get the correct leading nonanalytic (LNA) terms 
in the chiral expansion in any model where the pion-fermion coupling has a pseudoscalar 
(75) form [glHgl. 

The complete one-meson correction to the nucleon mass in this model can be calculated as 
the expectation value of A in the lowest-order wave function obtained from the homogeneous 
equation corresponding to Eq. (|l|). This can be written 

5M = ^A^, (16) 

where \1/ satisfies 

^ = GoK-^. (17) 

In models where the mesons are not introduced as elementary fields, some care should 
be taken with the quark-antiquark T- matrix, Eq. (|^). In the mesonic dressing of a single 
quark or one-loop corrections to meson properties, this T-matrix can be used as a meson 
propagator to generate a complete set of terms at next-to- leading order in 1/Nc where Nc is 
the number of colours [^,^,0]. These include terms which are linear in the coupling g^^, 
such as the exchange or Fock term corresponding to the Hartree self-energy of Eq. (H). For 
a diquark or a baryon, these linear terms are already included in the interactions between 
the quarks. Hence meson exchanges between two quarks which do not otherwise interact 
should be described by the subtracted T-matrix 

D^^ip) = D^^{p) - AigJ^^. (18) 

In the case of the pion, this has the form 

jj _ -83;n,(p-) 

l-2ig,n,(p^)- ^'"J 
This subtracted T-matrix should be used in the diagrams of Fig. ^, which form part of Fig. 

For a single quark (or a meson) the mesonic one-loop diagrams, including the Fock terms, 
form the leading corrections to the properties of these particles in a systematic expansion 
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in powers of l/Nc- In the case of baryons, things are not so clean. Both Hartree and Fock 
terms contribute to the baryon energy at leading order (A'^^). Meson exchanges between the 
quarks in a baryon are also of leading order, as are nonlinear terms arising from interactions 
among the mesons. In the large- A^^^ limit, this provides a justification for treating the baryon 
as a soliton, as in the Skyrme model. Such treatments have been quite widely applied to 
NJL-type models [p!2| , p!3| . They form a radically different approach to baryon structure from 
the one outlined here and it is still unclear which is the better for the case of interest. 



Nc = 3. Indeed a hybrid of the two has been proposed in Ref. [Q. Ultimately, numerical 
calculations of the size of the kernel based on Eq. (113) will indicate whether a perturbative 
treatment is adequate. 

If it is found that the shorter-range pieces of D are too strong to be treated perturbatively, 
then one possible way to improve the convergence of the approach may be to vary the 
strength of the subtraction term. This must be compensated by adjusting the strength of 
the interaction used in the bare equation for G. That would allow the strong short-range 
effects to be treated self-consistently in this equation, leaving the longer-range effects to act 
clS cl residual interaction. 

Note that the loop integral Eq. (p!0|) falls off at high momenta in this type of model. 
Consequently the subtracted T-matrix Eq. ( |I8| ) also falls off; the enhancements due to the 
Goldstone-boson nature of the pion only appear at low momenta. It is thus the underlying 
interaction responsible for spontaneous symmetry breaking which controls the short-range 
interaction between quarks, not the short-range part of Goldstone-boson-exchange as sug- 
gested by some authors . Provided that the quarks in a baryon can be assumed to interact 
in s-waves, only the strengths of the interactions in the scalar and axial diquark channels 
are relevant. 

The inclusion of interactions between the quarks while a meson is in flight. Fig. |^(b), 
is crucial to ensuring the correct energy denominators for the intermediate baryon-meson 
states. The binding energy of the nucleon depends on both the scalar and axial diquark 
couplings, whereas that of the A depends only on the axial one. The short-range interaction 
thus has an implicit spin-dependence, which splits the bare nucleon and A states. This means 
that, provided the scalar and axial diquark couplings are chosen to give the correct N-A 
splitting, the only intermediate states in the diagrams of Fig. H(b) with energy denominators 
of order m^r are pion-nucleon ones. These will generate the correct LNA terms in a chiral 
expansion, for example the order-m^ term in the nucleon mass 

A model of this type can therefore help to resolve some of the issues which have arisen 
in comparisons of the results of chiral quark models with those of the cloudy bag model 
In order to get the correct LNA terms (which are a signal that chiral symmetry 
has been correctly implemented) it is crucial that the full propagator G is used, allowing 
the quarks to interact while the meson is in flight. This means that the nucleon-plus-pion 
states necessary for the LNA terms are included, as in the cloudy bag model 0. However, 
in contrast to the elementary bag model (i.e. one defined strictly by the cloudy bag model 
Lagrangian [^), in our approach the spin-dependence of the interaction binding the "bare" 
baryons means that the A is not degenerate with the nucleon and so it does not contribute 
to the LNA terms in the nucleon mass. Although this problem is avoided in practical bag 
model calculations where the bare baryon in the baryon-plus-pion state is replaced by a 
physical one, such an approach is, apparently, equivalent to resumming a subset of mesonic 



8 



corrections to all orders in the Lagrangian-based description. Such a sum may give rise to 
conflicts with chiral symmetry at higher orders. In a related context, it has been noted that 
if ttTV loops are resummed in this way then one should also dress the vrA^ vertex with pion 
loops [0 . 



The chiral quark model of Glozman and coworkers [|6|, p7| , p0| includes pion loops which 
dress the quarks individually, effectively treating the quarks as non-interacting while the 
pion is in flight. In addition, retardation is neglected in the pion exchanges between quarks. 
These two approximations mean that the current treatments of this model do not generate 
the correct LNA terms. These would require the inclusion of interactions between the quarks 
while the meson is in flight, together with the corresponding retardation effects. 

Finally it is also worth noting that, in a model of the NJL type, a complete set of 
intermediate baryon states can be included in the meson-loop diagrams. This will make 
it possible to examine the role of excited baryons in the mesonic dressing. Although the 
contributions of individual excited states fall off with energy, there is an inflnite number of 
such states, which may be unbound as well as bound in a nonconfining model. The sum over 
these states would diverge in the absence of a cut-off. We therefore expect that they could 
make significant contributions to the mesonic dressing in the approach described here. In 
this way it will be possible to assess the accuracy of truncating to intermediate nucleon and 
A states only, as done in the original cloudy bag model and more recently in a model of 
the type discussed here p!6| , |T7 |. 

In this recent work, Oettel and coworkers |ll6| , p!7[| have calculated pion-loop contributions 
to nucleon properties within a Faddeev approach, with a covariant version of the approxi- 
mation used in the cloudy-bag model. In the original cloudy-bag model , pion loops were 
calculated keeping only vrA^ and vrA intermediate states and the baryons were treated as 
static. In Refs. |]TB|,|T^ these loops are calculated using fully covariant propagators for the 
nucleon and A but keeping only the poles corresponding to the lowest-energy solutions to the 
Faddeev equation in each channel. In esssence, this replaces their relativistic quark model 
by a model of elementary baryons and mesons. Oettel et al. also take care to include the 
tadpole diagrams required for the correct chiral behaviour. Their loop integrals are made 
finite by form factors at the meson-baryon vertices. These form factors can be regarded 
either as ad hoc regulators of the sort typically used in NJL models or as consequences of 
the internal structure of the hadrons. Although this approach respects chiral constraints 
and provides good estimates of the largest individual loop contributions, it does not give 
the complete loop correction since excited baryon states have been omitted. 



V. MESONIC DRESSING OF CURRENTS 

In order to calculate electromagnetic and weak properties of the nucleon, we need to 
construct the corresponding currents. At lowest order (no meson loops) the electromagnetic 
current operator to be used in a three-quark state can be obtained diagramatically by at- 
taching a photon everywhere in Eq. (0). The resulting expression for the matrix element of 
the current has the form 

J/^ = '^(f^ + K^')m. (20) 
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Here Fg denotes the sum of single-quark and diquark currents, and K^^ the interaction current 
where the photon is attached to the exchanged quark. This corresponds to the impulse 
approximation in the three-quark basis, properly subtracted to avoid double counting . 
The corresponding axial current can be found in the same way by making an axial-vector 
insertion instead of attaching a photon. Currents like Eq. pO| ) have been extensively studied 
in the NJL model as well as in QCD-motivated Faddeev |Q and SD approaches p4 



The same methods can also be used to calculate the matrix element of any operator which 
is bilinear in the quark fields, such as the pion-nucleon sigma commutator. 

At the next order, we need to consider the one-loop mesonic contributions to the currents, 
A^. These have a similar origin to the contributions to the kernel. A, and so can be obtained 
by attaching a photon everywhere in A, including the meson. The attachments to the 
composite mesons and diquarks are realised by direct couplings to quark currents in the qq 
and qq loops respectively. 

In addition we need to take into account the effects of mesonic fluctuations on the bound 
state wave function. If the full wave function is written as \& + 5^ , it satisfies the Bethe- 
Salpeter-type equation 



+ 5^) = Gq{K + A)(^ + 5m). 
The perturbation of the wave function can be expressed as p5 



5^ 



GbA 



2M 



dA ' 
\^)P^' vT/ 

9PM , 



|P|=A/ 



m. 



(21) 



(22) 



In the first term, Gb denotes the unperturbed Green function G with the bound-state pole 
subtracted off. This term incorporates configuration-mixing effects generated by meson- 
exchange between the quarks. The second term in this expression is just a wave-function 
renormalisation, analogous to that in the cloudy bag model. It involves the derivative of the 
meson-exchange part of the kernel with respect to the the total 4-momentum, evaluated at 
the unperturbed nucleon mass. 

The full one-loop contribution to the current matrix element is finally 



(23) 



Again, it is important that A^ contain the full propagator G, allowing the quarks to 
interact while the meson is in flight, see Fig. ^. As in the case of the nucleon mass, this 
generates the correct LNA terms in chiral expansions of the electromagnetic radii of the 
nucleon [|l|. In particular the isovector charge radius contains a term of order InrriTr and 
the isovector magnetic radius one of order m~^. In both cases these LNA terms contribute 
significantly to the full radii and so any calculation should include pion-nucleon loops with 
the correct energy denominators. 

Neglecting the interactions while the pion is in flight would lead to expressions for the 
currents like those in the quark model of Glozman and coworkers P, p0l , p6| , which do not 
contain the correct LNA terms. Truncating the sum over intermediate states to the nucleon 
and A only would lead to an approximation similar to the cloudy-bag model, with the correct 
LNA terms 
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FIG. 4. Mesonic corrections of three-quark (baryon) currents where the quarks interact while 
the meson is in flight. The dashed line denotes an external photon or other operator. 

VI. CONCLUSION 

We have presented a method for constructing a complete set of one-loop mesonic correc- 
tions to the nucleon in a chiral quark model, treated in a covariant manner. Our approach 
is based on the "gauging of equations" method which has previously been used to construct 
electromagnetic corrections to hadronic systems. It ensures that a complete set of one-loop 
corrections is generated, with no double counting. This completeness is essential if physical 
quantities calculated in this approach are to satisfy the constraints of chiral symmetry and 
gauge invariance. 

To illustrate the method, we have discussed in some detail its application to a particular 
covariant chiral quark model: the NJL model. The method starts from an unperturbed 
valence-quark Green's function, which in this case is obtained from a Faddeev equation. 
The complete set of mesonic insertions is used to generate the one-loop contribution to the 
kernel. The same method can also be used to construct photon-nucleon couplings in a way 
that respects current conservation. 

A central feature of this approach is the use of the unperturbed Green's function for three 
interacting valence quarks. This means that while a meson is in flight the quarks propagate 
in physical baryon states, as described at leading order in the meson-loop expansion. The 
fact that the NJL model gives rise to an N-A splitting at this order means that the only 
nonanalytic terms in the chiral expansion are generated by pion-nucleon states. This ensures 
that the method will give the correct leading nonanalytic terms in nucleon observables, in 
agreement with ChPT. This is particularly important for electromagnetic radii, where these 
nonanalytic terms make significant contributions. 

The usual treatment of the cloudy-bag model similarly gives the correct leading nonan- 
alytic terms in nucleon observables. However in models of the type discussed here, there 
is no need to truncate the unperturbed Green's function to just nucleon and A states. A 
complete set of baryon states can be summed over, regulated by the cut-off which forms an 
intrinsic part of any NJL-type model. 

In contrast, the chiral quark model of Glozman and coworkers effectively uses a free 
Green's function to describe the propagation of the quarks in meson-baryon intermediate 
states. It also neglects retardation effects on the exchange of pions between quarks. As a 
result this model, in the framework used so far, fails to give the correct nonanalytic terms in 
the chiral expansion. However this could be remedied by a treatment based on the approach 
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described here, or the corresponding time-ordered version. 
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